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ABSTRACT 

Boundary layer flow and turbulence transport 
analyses to study the Influence of the free-stream tur- 
bulence on the surface heat transfer rate and the skin 
friction around the stagnation point of a circular cyl- 
inder in a turbulent flow are presented. The analyses 
are formulated with the turbulent boundary layer equa- 
tions, the Reynolds stress transport equations and the 
k - c two-equation turbulence modeling. The analyses 
are used to calculate the time-averaged turbulence dou- 
ble correlations, the mean flow properties, the surface 
heat transfer rate and the skin friction with an Iso- 
tropic turbulence in the freestream. The analytical 
results are described and compared with the existing 
experimental measurements. Depending on the free- 
stream turbulence properties, the turbulence kinetic 
energy can Increase or decrease as the flow moves 
toward the surface. However, the turbulence kinetic 
energy Induces large Reynolds normal stresses at the 
boundary layer edge. The Reynolds normal stresses 
change the boundary layer profiles of the time-averaged 
double correlations of the velocity and temperature 
fluctuations, the surface heat transfer rate and the 
skin friction. The free-stream turbulence dissipation 
rate can affect the stagnation-point heat transfer rate 
but the Influence of the free-stream temperature fluc- 
tuation on the heat transfer rate Is Insignificant. 

INTRODUCTION 

The predictions of the stagnation-point flow prop- 
erties due to a circular cylinder In a turbulent flow 
are of practical Importance In engineering applications. 
In the Initial stage of the gas turbine blade design, 
the blade leading edge flow field Is usually simulated 
with the results of the stagnation-point flow studies. 
Estimation of the blade leading edge surface heat 
transfer Is an example. Much of the existing litera- 
ture Is focused on the study of the stagnation-point 
heat transfer rates. However, the Influence of the 
free-stream turbulence on the mean velocity, the mean 
temperature and the turbulence double correlations near 


the stagnation point Is critical to the prediction of 
the boundary layer flow development along the blade 
surface. Very limited results of the turbulence corre- 
lations and the mean flow properties within the stagna- 
tion-point flow were reported. This paper presents 
turbulence transport analyses of the velocity and the 
temperature fluctuations, the surface heat transfer 
rate and the skin friction within the stagnation-point 
flow of a circular cylinder in a turbulent freestream. 
The objective of this study Is to establish analyti- 
cally the stagnation-point boundary layer profiles of 
the mean velocity, the mean temperature and the time- 
averaged turbulence double correlations of the velocity 
and temperature fluctuations as functions of the turbu- 
lence Intensity, the turbulence length scale (turbu- 
lence dissipation rate) and the Reynolds number of the 
freestream. 

As In any turbulent flow analysis, modeling of the 
turbulence variation in the stagnation-point flow is a 
difficult task. The eddy viscosity concept was used In 
a turbulent boundary layer theory (Smith and Kuethe, 
1966) to analyze the surface heat transfer rate and the 
skin friction of the stagnation-point flow field. In 
order to Include the effect of the free-stream turbu- 
lence, the eddy viscosity was assumed to be propor- 
tional to the free-stream turbulence Intensity. 
Following the progress of the turbulence modeling tech- 
nique, Traci and Wilcox (1975) and Strahle (1985) ana- 
lyzed the boundary layer edge turbulence conditions. 
They used the two-equation turbulence models to 
describe the turbulence kinetic energy variation along 
the stagnation-point streamline. Hljlkata et al . 

(1982) added an equation for the anisotropic turbulence 
to the k - c turbulence modeling and numerically 
solved the Navler-Stokes equations for the Reynolds 
stress distributions along the stagnation-point stream- 
line. It Is generally agreed that a turbulence trans- 
port analysis using the Reynolds stress transport 
equations Is a better approach to study the Interaction 
between the turbulence and the mean flow strain. In a 
previous study (Wang and Yeh, 1987), Donaldson's turbu- 
lence modeling (Donaldson et al., 1972) and the k - e 
two-equation turbulence modeling (Launder and Spalding, 
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1974) were Implemented with the boundary layer flow 
equations to analyze the Influence of the fhee-stream 
turbulence on the surface heat transfer rate, the 
stagnation-point boundary layer profiles of the mean 
velocity, the mean temperature, and the time-averaged 
turbulence double correlations. The theoretical analy- 
ses of the previous study are used In this paper. How- 
ever, the existing analyses are Improved to model the 
variation of the time-averaged turbulence double corre- 
lations along the stagnation-point streamline and the 
boundary layer edge. The length scale for the turbu- 
lence closure assumptions are related analytically to 
the free-stream turbulence properties. The Improved 
analyses are also used to Investigate the Influence 
of the free-stream temperature fluctuation on the 
stagnation-point heat transfer rate. 

With the assumption of an Isotropic turbulence In 
the freestream, the analyses are used to calculate 
numerically the time-averaged turbulence double corre- 
lations and the mean flow properties within the bound- 
ary layer around the stagnation point. The surface 
heat transfer rate and the skin friction are also cal- 
culated. These computational results are compared with 
the measurements of existing experiments to validate 
the present analyses. 

NOMENCLATURE 

a mean velocity parameter 

cf skin friction factor 

q empirical constants, 1 - 1 to 4 

c u empirical constant 

D cylinder diameter 

E dimensionless turbulence dissipation rate 

F dimensionless mean flow velocity component 

G transformed velocity of V 

g turbulence longitudinal correlation 

H dimensionless mean temperature 


T mean temperature 

Tu free-stream turbulence Intensity 

t temperature fluctuation 

U.V mean velocity components 

u,v,w velocity fluctuations 

x.y.z physical coordinate system 

4 dimensionless length scale 

c turbulence dissipation rate 

n transformed coordinate along y direction 

0,9 dimensionless turbulence correlation 
A turbulence modeling length scale 

X turbulence longitudinal microlength scale 

p,v absolute and kinematic viscosity 

E transformed coordinate along x direction 

p density 

<f azimuth angle 

Subscripts 

e boundary layer edge condition 

F free-stream condition 

oe boundary layer edge condition at Initial station 
W surface condition 

0 condition at Initial station 

1 condition at location 1 

2 condition at location 2 of boundary layer edge 


K] dimensionless turbulence kinetic energy, 1 - 1 
to 3 

k turbulence kinetic energy 

a spatial distance 

Nu Nusselt number 

Pr laminar Prandtl number, 0.7 

R sum of dimensionless Reynolds normal stresses 

Re Reynolds number 

R x dimensionless Reynolds normal stress along x 

direction 

R y dimensionless Reynolds normal stress along y 
direction 

R z dimensionless Reynolds normal stress along z 

direction 


S,t dimensionless Reynolds shear stress 


Superscript 

time-averaged quantity 
THEORETICAL FORMULATION 


The theoretical model of the flow field Is shown 
schematically In Fig. 1. The stagnation point 0 Is 
taken as the coordinate-system origin with x and y 
coordinates parallel to and normal to the cylinder sur- 
face, respectively and z coordinate along the cylin- 
der axis. The diameter of the cylinder Is D and Its 
axis Is perpendicular to the free-stream mean velocity 
Vp. The mean flow Is assumed to be steady, Incompress- 
ible, and two-dimensional. The mean velocity compo- 
nents are U and V along the x and y directions. 
The turbulent velocity fluctuations along the x, y, 
and z directions are denoted by u, v, and w. The 
free-stream Reynolds number Is defined as Re - VpD/v, 
where v Is the kinematic viscosity. The free-stream 
turbulence Is characterized by Its turbulence Intensity 


Tu 
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c F . The free-stream mean temperature Is Tp. The cyl- 
inder surface Is at a constant temperature Ty. 

The following hypotheses are used In developing 
the theoretical formulation. A boundary layer flow Is 
assumed to exist along the surface around the stag- 
nation point. The mean velocity In the external flow 
Is represented by the result of the Invlscld flow anal- 
ysis around a circular cylinder. The turbulence Is 
distorted as the flow moves toward the stagnation point 
and the turbulence conditions, different from those In 
the freestream, are Imposed along the boundary layer 
edge. The turbulence boundary layer edge conditions 
induce turbulence momentum and thermal fluxes within 
the boundary layer. These turbulence fluxes Influence 
the boundary layer mean flow properties, the surface 
heat transfer rate and the skin friction. Based on the 
above hypotheses, the following theoretical models are 
proposed to Investigate the momentum and thermal Inter- 
action within the stagnation-point flow field. 

Boundary Laver Flow Analysis 

Similar to the analyses In the previous study 
(Wang and Yeh, 1987), the mean flow properties within 
the boundary layer are described by the two-dimensional 
boundary layer equations. The turbulence flux terms In 
momentum and enthalpy equations are modeled with the 
time-averaged turbulence transport equations. The mean 
flow properties and the turbulence double correlations 
are nondlmenslonal lzed and they are defined as 



The physical coordinates x and y are also 
transformed Into the dimensionless coordinate system !• 
and n by using the following relations: 
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and n 
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Momentum equation, 
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It Is further assumed that the turbulence does 
not Influence the mean velocity within the external 
flow region and Its mean velocity Is represented with 
the two-dimensional Invlscld flow solution due to a 
circular cylinder In a uniform freestream. The mean 
velocity components are 
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Therefore, the governing equations for the mean 
flow and the turbulence become 

Continuity equation, 

2E + |^ + F - 0 (3) 
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where R » R x + R y + R z - 
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The above governing equations were derived In the 
previous work (Wang and Yeh, 1987) by following the 
theoretical analyses by Donaldson et al . (1972) and 
Beckwith and Bushnell (1966). A length scale A was 
used In the closure assumptions for the time-averaged 
higher-order turbulence correlations. In this study, 
the following linear relations are assumed for the 
length scale within the boundary layer, 

A = y, for 0 < y < y p 

( 11 ) 

A ■ \ - v(v^l) f ° r 


where y p Is a very small distance (» 0.02 y e ) from 
the surface. The subscript p denotes the pivot-point 
locations . 

Along the stagnation-point streamline, the upstream 
limit Hi of the external flow region Is defined by 
matching the mean flow velocity, I.e., Vj - V p at 

n » nj- The downstream limit n 2 defined by matching 

the turbulence conditions at the boundary layer edge. 

The boundary layer edge conditions of the turbu- 
lence double correlations are analyzed with the turbu- 
lence transport equations. For this purpose, the 
fol lowl ng conditions 
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Where K] , < 2 , K 3 , x, 6 , and A are the dimensionless 
turbulence double correlations and length scale. They 
are defined as 
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-uv 
av ’ 



are assumed at the locations near the boundary layer 
edge . 

Nonsllp conditions 
T 

F = 0, H . f. , R x . 0. R y - 0, R z - 0. 


with K - K ] + K 2 + K 3 . 


The free-stream turbulence properties are the 
turbulence conditions at ni • Therefore, at . 



S ■ 0, and 0 ■ 0 (13) 


e - e F (19) 


are Imposed along the cylinder surface n • 0. Based 
on the above boundary conditions, Eqs. (3) to (10) are 
used In this study to model the mean flow properties 
and the turbulence double correlations. 


Stagnation-Point Streamline Turbulence 

The turbulence double correlations along the 
stagnation-point streamline within the external flow 
region are modeled with Donaldson's turbulence model- 
ing technique. With the mean flow velocity assumption, 
Eq. (2), the turbulence transport equation becomes 
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The turbulence conditions at n 2 were described 

In Eq. (11). However, the Reynolds normal stress 
gradients at n 2 are replaced by their corresponding 
Reynolds normal^stresses . The following analysis Is 
used to determine these normal stress components. 

Applying the assumptions In Eq. (12) to the 
turbulence transport Eqs. ( 6 ) to (10) results In the 
following expressions for the gradients of the tur- 
bulence correlations along the E direction 
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It Is assumed that, within the boundary layer 
thickness In the vicinity of the stagnation point, the 
Reynolds normal stresses are Independent of v 
Therefore, 
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are found at v * 0 . 

Substituting the above relations Into Eqs. (20) 
to ( 22 ) and rearranging the equations leads to the 
fol lowing expressions 
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Equations (29) to (31), dt/dq - 0, and de/dn = 0 
are the final forms of the turbulence boundary condi- 
tions at n 2 - With the turbulence boundary conditions 
at q, and^ n 2 . Ec is. (14) to (18) are used to model 
the turbulence^double correlations at locations between 
n 1 and n 2 - A constant A oe Is also assumed In this 

region. This length scale Is determined by the require- 
ment that 
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The k - c two-equation turbulence modeling Is 
used to define K oe . Based on the previous analysis 
(Wang and Yeh, 1987), the k and e equations are 
written as 
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with 
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for the Reynolds normal stresses at v 0 ° 0 and 
n » n 2 ■ To be conslstant with the definitions of the 

variables In Eqs. (14) to (18), the above relations are 
rewritten as 
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The boundary conditions are 
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at n. 


and 

dK _ 

^ - 0 at n 2 . 

The empirical constants are ci » 1.44, C 2 - 1.92, 

C 3 - 1.00, C 4 - 1.30, and cp - 0.09. 

The above boundary value problem defines the 
turbulence kinetic energy K in n 1 > n 2 n 2 along 

the stagnation-point streamline. The previous work 
(Wang and Yeh, 1987) also Indicated that the boundary 
layer edge was located at n 2 ~ 5. Once the turbu- 
lence kinetic energy Ko e Is obtained, K] >oe , K2 i 0 e> 
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K 3 oe , and A oe can be found by solving Eqs. (29) 
to’ (32) . 

The previous analysis (Wang and Yeh, 1987) Is also 
used to model the Initial profiles of the turbulence in 
n ? > n > 0. The turbulence correlation gradient terms 
in ^-direction are also assumed to be Independent of 
n and they are represented by Eqs. (20) to (24). Oue 
to symmetry In the mean flow, 3F/3E - 0 and 3H/3E - 0 
are true along the stagnation-point streamline. Sub- 
stituting these mean flow streamwlse gradients Into 
Eqs. (3) to (10) results In a set of ordinary differ- 
ential equations which describe the variations of the 
mean flow properties and the turbulence double corre- 
lations in n 2 > n > 0. Equations (12) and (13) 

express the boundary conditions at n » 0 and n 2 - 

Boundary Layer Edge Turbulence 

Equations (3) to (10) are also used to model the 
turbulence double correlations within the boundary 
layer flow around the stagnation point. The Initial 
profiles of the mean flow properties and the turbu- 
lence double correlations within the boundary layer 
thickness along the stagnation-point streamline were 
defined from the analysis In the previous sections. 

The surface boundary conditions were also described by 
Eq. (13). The boundary layer edge mean flow properties 
are F - 1 and H - 1. Thus, the boundary layer edge 
turbulence conditions are required to form a complete 
set of the mean flow and turbulence boundary conditions 
for Eqs. (3) to (10). Analysis of the boundary layer 
edge turbulence conditions Is described In the 
following. 

The k - c two-equation turbulence modeling Is 
used again to predict the turbulence kinetic energy 
along the boundary layer edge. Following the previous 
analysis (Wang and Yeh, 1987), the k and e equa- 
tions are written as 


(Wang and Yeh, 1987). These existing numerical compu- 
tational procedures were performed to verify the 
present theoretical formulation. This verification Is 
also based on existing experiments. Several sets of 
the experimental data on the effect of the free-stream 
turbulence on the heat transfer rate from a heated 
cylinder placed normal to alrstream were reported by 
Smith and Kuethe (1966), Lowery and Vachon (1975), and 
Zukauskas and Zlugzda (1985). Honeycomb structure, 
damping screens or turbulence Inducing grids at various 
locations upstream of the test cylinder provided varia- 
tions In the free-stream turbulence conditions. The 
measurements of the two-point turbulence longitudinal 
double correlations In the free stream were also 
reported by Lowery and Vachon (1975). The two-point 
turbulence correlations provide information regarding 
the free-stream turbulence dissipation rate. There- 
fore, their test conditions and results are first used 
to validate the present theoretical formulation. 
Assuming the free-stream turbulence Is Isotropic, the 
free-stream turbulence kinetic energy and turbulence 
dissipation rate are described (Hlnze, 1975) by 

3v? 30w? 

kp = &nd Cp ■ — 2 — (37) 

Xp 

The longitudinal two-point turbulence correlation 
g, the spatial distance 2 ., and the turbulence longi- 
tudinal microlength scale Xp are related by the 
following relations for small 1 
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(35) 

(36) 


The analysis along the stagnation-point streamline 
provides the K and E Initial conditions at (E 2 .n 2 )- 

Equations (35) and (36) are solved by numerical Itera- 
tion for the boundary layer edge turbulence kinetic 
energy within the small region <p < 10 °. Equations (20) 
to (24) are also Integrated numerically to calculate the 
downstream boundary layer edge turbulence kinetic energy 
(or the corresponding Reynolds normal stresses). The 
length scale A e Is not necessarily constant along the 
boundary layer edge. Through numerical Iteration, this 
local length scale A e Is found such that the sum of 
the Reynolds normal stresses converges to the local 
turbulence kinetic energy obtained from eqs. (35) 
and (36). The local Reynolds normal stresses are the 
boundary layer edge turbulence conditions. 


NUMERICAL COMPUTATION 


Most of the numerical computational procedures 
used In this study were described In the previous work 


The turbulence longitudinal microlength scales are 
obtained by curve-fitting the experimental data with 
Eq. (38). Based on the free-stream mean velocity, the 
free-stream turbulence Intensity, and the turbulence 
longitudinal microlength scale, the present analysis 
Is used to calculate numerically the turbulence double 
correlations and the mean flow properties within the 
stagnation-point flow field. The surface heat trans- 
fer rate and the skin friction are then determined 
from the mean velocity and the mean temperature 
profiles. 

Summary of the computational sequences associated 
with the present theoretical formulation are presented 
here. The computational sequences are: 

(1) The k - c two-equation turbulence modeling, 
Eqs. (33) and (34), were solved for the turbulence 
kinetic energy K In n 1 > n > n 2 - The turbulence 

kinetic energy at n 2 was substituted Into Eqs. (29) 
to (31) to solve for^the turbulence kinetic energy 
K, antl by numerical Iteration In 

I , 06 c. 1 06 j ) 06 

terms of the length scale A 0 p. The length scale Ao e 
was chosen such that the resulting turbulence double 
correlations satisfied Eq. (32). 

(2) Turbulence kinetic energy In n 1 > n > n 2 and 

the length scale A oe were used to linearize the finite 
difference equations corresponding to Eqs. (14) to (16). 
The finite difference equations were obtained by using 
central finite difference approximation. Since the 
Reynolds normal stresses were specified at ru 2 n > n 2> 
these finite difference equations were solved with the^ 
method of successive substitution. These Reynolds nor- 
mal stress distributions were used to calculate z and 
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0 in n 1 2 n 2 n 2 from the finite difference equa- 
tions corresponding to Eqs. (17) and (18). The free- 
stream condition x - 0 and an arbitrary boundary 
condition 9 « 9p were imposed at n 1 In order to 

study the effect of the correlation of the free-stream 
velocity and temperature fluctuations on the surface 
heat transfer rate. These finite difference equations 
were solved by using the method of the matrix Inver- 
sion. Variable grid dimension was used In numerical 
computation. Two hundred grid points were used. The 
grid was tightened near ni and q. to ensure the 
boundary conditions. ' '' 

(3) The turbulence double correlations at (E 2 .n 2 * 

were substituted Into Eqs. (20) to (24) to determine 
the E-dlrectlon gradients. The mean flow conditions 
were calculated at <? - 0.01° location. Introducing 
the E-dlrection gradients into Eqs. (3) to (10) led 
to a set of ordinary differential equations with the 
boundary conditions specified at n - 0 and n 2 - 

These equations were then solved numerically for the 
mean flow properties and the turbulence correlations. 

It was also found that the choice of <p( » 0*) other 
than 0 . 01 ° did not significantly change the results 
of the numerical computations. This computational 
step and step 2 calculated the turbulent flow proper- 
ties along the stagnation-point streamline. 

(4) Based on the Reynolds normal stresses at 
(E 2 ,n 2 ). Eqs. (35) and (36) were solved numerically 

for the turbulence kinetic energy along the boundary 
layer edge. A fourth-order Runge-Kutta numerical 
method was used to Integrate Eqs. (20) to (24) for the 
downstream turbulence double correlations. Very small 
step size In ( 2 ? corresponding to A? - 0 . 001 °) was 
used. The numerical Integration at any location E 
was performed Iteratively In terms of the length scale 
A e . The Iteration procedure was terminated when the 
sum of the Reynolds normal stresses converged to the 
turbulence kinetic energy previously obtained from 
Eqs. (35) and (36). Since x » 0. and 6 - 0. were 
calculated at (E 2 ,n 2 ). Eqs. (23) and (24) showed that 
x - 0 . and 9 - 0 . also occurred along the downstream 
boundary layer edge. 

(5) Finally, the properties of the mean flow and 
the turbulence double correlations within the boundary 
layer flow around the stagnation point (<p S 10 °) were 
calculated from Eqs. (3) to (10) by using an existing 
downstream marching numerical computational scheme. The 
Nusselt number Nu and the skin friction factor Cf 
were then calculated from the mean velocity and the 
mean temperature profiles. They are defined as 



The mean velocity and mean temperature gradients at 
the surface were calculated numerically from the mean 
velocity and mean temperature profiles which were previ- 
ously obtained from the boundary layer flow solutions. 

RESULTS AND DISCUSSIONS 

The numerical computational procedures were per- 
formed for several cases with different free-stream 
turbulence conditions. Examples of the freestream tur- 
bulence longitudinal double correlations and the turbu- 


lence microlength scale are shown In Fig. 2. Some of 
the present computational results are described In the 
following sections and these results are also compared 
with available experimental measurements. 

Turbul ence 

The variations of the dimensionless Reynolds nor- 
mal stresses along the stagnation-point streamline are 
shown In Figs. 3(a) and (b). Rapid changes In the 
normal stresses occur near ni ■ for most cases, the 
normal stress components Ki and K 3 decrease and <2 
Increases along the stagnation-point streamline within 
the external flow region. Within the boundary layer, 
the Reynolds normal stresses decrease continuously In 
the direction toward the surface. Large K 2 and small 
K] and <3 were predicted at locations close to the 
surface. Figure 3(b) also shows that, with Tu = 0.028, 
Re - 2.52x10 s , and Xp/D =. 0.0079, the normal stresses 
amplify along the stagnation-point streamline within the 
external flow and the maximum normal stresses occur at 
the boundary layer edge. 

_ The profiles of the dimensionless shear stress 
-uv/va and the dimensionless mean velocity U/U e 
within the boundary layer thickness along the 
stagnation-point streamline are shown In Figs. 4(a) 
and (b). With the prescribed free-stream turbulence 
longitudinal microlength scales, these results show 
that large shear stress was predicted within the 
boundary layer for the case with large Tu or Re. 
Increasing Tu or Re Increases U/U e within the 
region close to the surface and decreases U/U e at 
locations near the boundary layer edge. The existing 
measurements (Kestln and Richardson, 1969) of U/U e 
with Tu - 0.05 and Re * 2.30x10 s were also plotted 
In Fig. 4(a) to compare with the present analytical 
results of U/U e . Following the test conditions of 
Lowery and Vachon (1975), Re - 1.11x10 s and 
Xp/D - 0.0086 were used In the present numerical 
calculation and the present analysis predicted a 
smaller U/U e In the region near the surface and a 
larger U/U e at locations near the boundary layer 

edge. — 0 5 

The dimensionless turbulence correlation -vt/kp Tp 

and mean temperature T/Tp along the stagnation-point 

streamline are shown 1 n_F 1 gs . 5(a) and (b). The tur- 
bulence correlation _^vt decays rapidly within the 
external flow and -vt - 0 occurs at the boundary 
layer edge. This turbulence correlation appears again 
within the boundary layer. It Increases to a maximum 
value near the surface and then decreases to the non- 
slip surface condition. With the prescribed turbulence 
longitudinal microlength scale, 1 ncreas 1 ng_the free- 
stream turbulence Intensity changes the -vt and T 
profiles within the boundary layer and Increases the 
mean temperature gradient at the surface. Therefore, 
the turbulent freestream Induces a larger surface heat 
transfer rate than that with a laminar free-stream 
condition. 

Analytical results of the Reynolds normal stresses 
and the length scales for the turbulence closure 
assumptions along the boundary layer edge are shown in 
Figs. 6 (a) and (b). The Reynolds normal stresses 
remain constant near the stagnation point and then 
Increase continuously at the downstream locations. 

The turbulence Is anisotropic and the component K 2 t e 
Is the largest among the Reynolds normal stresses. 
Constant length scale A e was predicted for the cases 
with large Tu or large Re. The length scale 
Increases slightly at the downstream locations for the 
case with a small Tu. In the analysis of the boundary 
layer edge turbulence, the turbulence kinetic energy 
was determined by solving Eqs. (33) and (34). The 
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diffusion terms are retained In these equations. Exam- 
inations of the present computational results indi- 
cate that the production and the diffusion terms In 
Eqs . (33) and (34) are the dominant terms In the turbu- 
lence development. The diffusion term is very small 
except at the locations very near the stagnation point. 
Therefore, the present computational results from 
Eqs. (33) and (34) are consistent with the boundary 
layer flow assumption. 

The computational results of the Reynolds normal 
stresses within the boundary layer at <f> - 0°. 5°, and 
10° are shown In Figs. 7(a) and (b). These results 
indicate that large Reynolds normal stresses occur at 
the boundary layer edge and the normal stresses 
decrease continuously along the direction toward the 
surface. The components K] and Kn Increase as the 
flow moves downstream. Within a small region near the 
surface, the <2 component decreases at the locations 
near the stagnation point and It Increases again at 
the downstream locations. These computational results 
of the normal stress profiles are different from the 
existing measurements of the normal stress profiles 
within a flat plate turbulent boundary layer 
(Klebanoff, 1955). 

The Reynolds shear stress profiles within the 
boundary layer at - 0*. 5°, and 10° were plotted In 
Figs. 8 (a) and (b). These profiles show_ that the shear 
stress Increases from the zero value at 'the boundary 
layer edge to a maximum value at n - 1 and then 
decreases to the boundary value at the surface. The 
shear stress also Increases significantly within the 
boundary layer at the downstream locations. 

The analytical results of the profiles of the 

dimensionless time-averaged correlation -vt/kp'^Tp 

within the boundary layer at <p - 0°, 5°, and 10* are 
shown In Figs. 9(a) and (b). These profiles are simi- 
lar to the shear stress profiles In Figs. 8 (a) and (b). 
The turbulence correlation also Increases as the flow 
moves downstream. With a large Tu ln_the free stream 
(Fig. 9(a)), small variations In the -vt profiles 
occur within the downstream boundary layer. However, 
the profiles In Fig. 9(b) show that, with a small Tu 
in the free stream, the -vt correlation Increases 
within the boundary layer at the downstream locations. 

2 

Since k = Tu , these profiles also Indicate that a 
large Tu or a large Re In the free stream can 
induce a large level of the -vt correlation within 
the boundary layer. 

Surface Heat Transfer and Skin Friction 

The present analytical results of the Nusselt 
number and the skin friction factor within 10 * of the 
stagnation point are shown In Figs. 10(a), and (b). 

These results are compared with the stagnation-point 
Nusselt number from the existing experiment (Lowery 
and Vachon, 1975) and the skin friction factor from 
the laminar boundary layer flow analysis (Schllchtlng, 
1968). The present analytical results of the 
stagnation-point Nusselt number are In good agreement 
with the experiment, with Tu < 0.055. For the case 
of Tu - 0.110, the present analyses predict a lower 
level of the stagnation-point Nusselt number than the 
measurement. Since similar level of the Nusselt number 
was also measured by O'Brian and Vanfossen (1985), dif- 
ferent length scale may be required In the analyses to 
calculate these experimental stagnation-point Nusselt 
numbers at Tu - 0.110. The present analytical results 
also show that the Nusselt number Increases slightly at 
locations near the stagnation point and then remains 
constant at downstream locations. However, the skin 
friction factor Increases continuously at the down- 
stream locations. These surface phenomena are consist- 


ent with the variations of the boundary layer turbulence 
profiles (Figs. 4 and 52- Although the turbulence cor- 
relations -uv and -vt have decreased to zero at the 
boundary layer edge, Reynolds normal stresses occur 
within the boundary layer. These normal stresses Inter- 
act with the mean velocity and thejnean temperature to 
produce turbulence correlations -uv and -vt within 
the boundary layer. Therefore, depending on the free- 
stream velocity turbulence, the present analysis pre- 
dicts larger surface heat transfer rate and skin 
friction than those with laminar flow conditions. 

To show the effect of the free-stream turbulence 
longitudinal microlength scale on the surface heat 
transfer rate, some computational results for the 
cases with Re - 2 . 52 x 105 , Tu - 0.028, and different 
Xp/D - 0.0053 and 0.0079 are compared In Figs. 11(a) 
to (d). With a small Xp/D (« 0.0053) In the free- 
stream turbulence, small amplification In <2 ancl 
attenuations in Kj and K 3 were predicted along the 
stagnation-point streamline (Fig. 11(a)). With a large 
Xp/D (» 0.0079), the Reynolds normal stresses amplified 
along the stagnation-point streamline In the external 
flow and these_normal stresses Induced large turbulence 
correlation -vt within the boundary layer and mean 
temperature gradient at the surface (Fig. 11(b)). 

The variations of the turbulence kinetic energy 
along the stagnation-point streamline are also shown 
In Fig. 11(a). Depending on the free-stream turbulence 
microlength scale, the turbulence kinetic energy can 
Increase or decrease as the flow moves toward the stag- 
nation point. With a small Xp/D, the free-stream tur- 
bulence has a large turbulence dissipation rate. This 
large dissipation rate Is the dominant factor In the 
development of the turbulence and the turbulence kinetic 
energy Is reduced. A smaller K e , as It Is compared 
with Kp, occurs along the boundary layer edge and this 
small turbulence kinetic energy Induces a low level of 
the turbulence correlation -vt within the boundary 
layer. Therefore, small mean temperature gradient and 
heat transfer rate occur at the surface (Figs. 11(b) 
and (c)). For the case with a large free-stream turbu- 
lence microlength scale, the production term Is the 
dominant factor In the turbulence development and the 
turbulence kinetic energy Is Increased. Therefore, a 
large turbulence kinetic energy occurs at the boundary 
layer edge and this turbulence kinetic energy_ 1 _nduces 
a high level of the turbulence correlation -vt within 
the boundary layer. Therefore, large mean temperature 
gradient and heat transfer rate occur at the surface 
(Figs. 11(b) and (c)). With Xp/D = 0.0079 In the 
free-stream turbulence, the present analysis predicted 
the stagnation-point Froessllng number reported by 
Lowery and Vachon (1975). 

Smith and Kuethe (1966) and Zukauskas and Zlugzda 
(1985) also performed the experiments of the stagnation- 
point surface heat transfer rates due to circular 
cylinders In turbulent freestream. These existing 
experimental results have been used by many researchers 
to verify their analyses of the stagnation-point flow. 
The present analysis was also used to predict these 
existing measurements of the stagnation-point heat 
transfer rates. Since details of the free-stream tur- 
bulence length scales were not reported In these 
studies, the present analysis could not be used to 
calculate the turbulence kinetic energy along the 
stagnation-point streamline within the external flow. 
However, overall heat transfer results could be 
attempted. The assumption that the boundary layer edge 
turbulence kinetic energy Is proportional to the free- 
stream turbulence kinetic energy was made to perform 
the boundary layer flow analysis. The computational 
results, based on three different K e /Kp ratios, of 
the stagnation-point heat transfer rates are compared 
with the various measurements In Fig. 11(d). These 
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results are presented In terras of the correlation 
parameters TuRe^'S and NuRe - ^-^. With K e /Kp - 1, 
the present analytical results agree with the experi- 
ments and the surface heat transfer correlation 
reported by Smith and Kuethe (1966). A smaller K e , as 
it Is compared with Kp, was required In the present 
analysis to predict the existing measurements of Lowery 
and Vachon (1975) and Zukauskas and Zlugzda (1985), 
especially for the cases with large TuReO-5. With a 
high turbulence Intensity In the freestream, the turbu- 
lence Is most likely to be anisotropic and Is dominated 
by the turbulence along the Vp direction. The free- 
stream turbulence has a lower turbulence kinetic energy 
than that of an Isotropic turbulence and 1 t_has a large 

2 

turbulence dissipation rate due to large Vp. With a 

large free-stream turbulence dissipation rate, the tur- 
bulence kinetic energy Is reduced as the flow moves 
along the stagnation-point streamline. A small turbu- 
lence kinetic energy occurs at the boundary layer edge 
and this turbulence kinetic energy Induces a low sur- 
face heat transfer rate. Therefore, K e < Kp Is 
required In the present boundary layer flow analysis to 
predict the experiment with a large value of TuRe^-5 
in the freestream. It Is also Interesting to see that 
the present analysis can predict the experimental 
stagnation-point Froessllng number (Zukauskas and 
Zlugzda, 1985) which Is out of the range of the exist- 
ing correlations by Smith and Kuethe (1966) and Lowery 
and Vachon (1975). This result Is significant since 
the experiment by Zukauskas and Zlugzda (1985) was 
pressurized. The freestream Reynolds number was large 
and the boundary layer thickness over the cylinder sur- 
face was very thin In their tests. 

The analysis of the length scale for the turbulence 
closure assumptions Is another aspect of the turbulence 
modeling. In general, two different length scales are 
required. The length scales are related to each other 
through additional empirical constants and a turbulence 
Reynolds number (Donaldson et al., 1972). This type of 
analysis on the length scale was not pursued In this 
study. Based on the discussion of the length scales In 
the work by Donaldson et al . (1972), simple relations, 
Eq. (11), for a single length scale was used In the 
present analysis. However, the effect of the free- 
stream turbulence on the length scale was considered 
through the analysis of the boundary layer edge turbu- 
lence conditions. Equation (11) enabled the present 
analysis to predict satisfactorily the existing meas- 
urements of the stagnation-point heat transferrates. 

CONCLUSIONS 

Boundary layer flow and turbulence transport 
analyses were formulated to study the development of 
the time-averaged double correlations of the velocity 
and temperature fluctuations and the mean flow pro- 
perties In the stagnation-point flow of a circular 
cylinder In a turbulent freestream. Steady and 
Incompressible turbulent boundary layer flow equa- 
tions, the k - c two-equation turbulence modeling 
and the Reynolds stress transport equations were used 
to formulate the theoretical models. The boundary 
layer edge turbulence conditions and the length scale 
for the turbulence closure assumptions were also 
obtained analytically by matching the turbulence 
boundary conditions. With the assumption of an Iso- 
tropic free-stream turbulence, the analyses were used 
to calculate numerically the surface heat transfer 
rate, the skin friction and the boundary layer profiles 
of the mean velocity and the mean temperature In the 
stagnation-point flow. The present analytical results 
were then compared with the existing experiments to 
validate the theoretical models In the analysis. 


The results of this study indicate that: 

1. The time-averaged double correlation of the 
velocity and temperature fluctuations decays rapidly 
along the stagnation-point streamline within the 
external flow region. The turbulence correlation of 
the velocity and temperature fluctuations Is very small 
at the boundary layer edge. Therefore, the free-stream 
temperature fluctuation does not appear to penetrate 
Into the boundary layer and has no effect on the sur- 
face heat transfer rate. 

2. With turbulence In the free stream, the 
Reynolds normal stresses may Increase or decrease as 
the flow moves toward the surface along the stagnation- 
point streamline. However, large Reynolds normal 
stresses usually occur at the boundary layer edge. 

These Reynolds normal stresses Induce the shear stress 
-uv and the turbulence thermal flux -vt within the 
boundary layer. These turbulence fluxes change the 
boundary layer mean velocity and mean temperature pro- 
files and, thus, the surface heat transfer rate and the 
skin friction. Including the effect of the free-stream 
turbulence, the present analysis predicts larger levels 
of the surface heat transfer rate and the skin friction 
than those with laminar free-stream conditions. 

3. The free-stream turbulence dissipation rate can 
Influence the stagnation-point heat transfer rate. 

When the free-stream turbulence dissipation rate Is 
large, the turbulence kinetic energy Is reduced as the 
flow moves toward the boundary layer edge. Due to the 
small turbulence kinetic energy at the boundary layer 
edge, a small surface heat transfer rate occurs beneath 
the stagnation-point flow. The free-stream turbulence 
dissipation rate Is found to be a factor In the differ- 
ence of the stagnation-point heat transfer rates 
reported In the existing studies. 

4. With the prescribed Isotropic free-stream tur- 
bulence conditions, the present analysis calculates the 
existing experimental measurements of the stagnation- 
point heat transfer rate when the free-stream Reynolds 
number Is In the order of 10^ and the free-stream tur- 
bulence Intensity Is <10 percent. The present analysis 
can be used to calculate the turbulence and mean flow 
properties when the free-stream turbulence Is aniso- 
tropic. However, It requires the free-stream Reynolds 
normal stresses and the free-stream turbulence dissipa- 
tion rate. 

REFERENCES 

Beckwith, I.E., and Bushnell, D.M., 1968, "Detailed 
Description and Results of a Method for Computing Mean 
and Fluctuating Quantities In Turbulent Boundary 
Layers," NASA TN D-4815. 

Donaldson, C.D., Sullivan, R.D., and Rosenbaum, H., 
1972, "A Theoretical Study of the Generation of 
Atmospheric-Clear Air Turbulence," AIAA Journal , 

Vol. 10, No. 2, pp. 162-170. 

HI j 1 kata. K., Yoshlda, H., and Mori, Y., 1982, 
"Theoretical and Experimental Study of Turbulence 
Effects on Heat Transfer Around the Stagnation Point 
of a Cylinder," Heat Transfer 1982 , Vol. 3, U. Grlgull, 
et al, eds., Hemisphere Publishing Corp., Washington, 
D.C., pp. 165-170. 

Hlnze, J.O., 1975, Turbulence . 2nd edition, McGraw- 
Hill, New York. 

Kestln, J., and Richardson, P.D., 1969, "The Effects 
of Free-Stream Turbulence and of Sound Upon Heat 
Transfer," ARL-69-0062, Wrlght-Patterson AFB, Dayton, 
OH. (Avail. NTIS, AD-691497). 

Klebanoff, P.S., 1955, "Characteristics of Turbulence 
In a Boundary Layer with Zero Pressure Gradient," NACA 
Rept. 1247. 


9 














TURBULENCE LONGITUDINAL CORRELATION 



FIGURE 2. - FREE-STREAM TURBULENCE LONGITUDINAL MICRO- 
LENGTH SCALE/ REYNOLDS NUMBER, Re = 2.52 X 10 5 . 
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(B) Re = 2.52 x 10 5 . 

FIGURE A. - VARIATIONS OF MEAN VELOCITY AND REYNOLDS SHEAR 
STRESS WITHIN BOUNDARY LAYER ALONG THE STAGNATION-POINT 
STREAMLINE. 
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MEAN TEMPERATURE, T/T p 

(B) Re = 2.52 x 10 5 . 

FIGURE 5. - VARIATIONS OF MEAN TEMPERATURE AND THE TUR- 
BULENCE CORRELATION WITHIN BOUNDARY LAYER ALONG STAG- 
NATION-POINT STREAMLINE. 
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LENGTH SCALE 




(A) Re = 1.11 X 10 5 . 
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(B) Re = 2.52 X 10 5< 


FIGURE 6. - VARIATIONS OF REYNOLDS NORMAL STRESSES AND LENGTH SCALE FOR 
TURBULENCE CLOSURE ASSUMPTION ALONG BOUNDARY LAYER EDGE. 
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DISTANCE NORMAL TO SURFACE 




(B) Re = 2.52 X 10 5< 

FIGURE 7. - REYNOLDS NORMAL STRESS PROFILES WITHIN BOUND- 
ARY LAYER AT DOWNSTREAM LOCATIONS. 
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FIGURE 8. - REYNOLDS SHEAR STRESS PROFILES WITHIN BOUND- 
ARY LAYER AT DOWNSTREAM LOCATIONS. 
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(B) Re = 2.52 x 10 5 . 

FIGURE 9. - PROFILES OF TURBULENCE CORRELATION - vt 
WITHIN BOUNDARY LAYER AT DOWNSTREAM LOCATIONS. 
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TURBULENCE KINETIC ENERGY. K 
(A) REYNOLDS NORMAL STRESSES AND TURBULENCE KINETIC 
ENERGY. 
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(B) MEAN TEMPERATURE AND TURBULENCE CORRELATION. 
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(C) BOUNDARY LAYER EDGE TURBULENCE KINETIC ENERGY AND (D) STAGNATION-POINT HEAT TRANSFER CORRELATIONS. 

FROESSLING NUMBER AT LOCATIONS NEAR STAGNATION POINT. 


FIGURE 11. - EFFECT OF FREE-STREAM TURBULENCE LENGTH SCALE ON TURBULENCE CORRELATIONS, MEAN TEMPERATURE PROFILES, AND STAGNATION- 
POINT HEAT TRANSFER CORRELATIONS. 


> 


20 



NASA 

National Aeronautics and 
Spaca Admlnlatratlon 

1. Report No. 


Report Documentation Page 

I 2. Government Accession No. I 3. Recipient's Catalog No. 


NASA TM-1 


iMl 


4. Title and Subtitle 

Correlations of Velocity and Temperature 
Fluctuations In the Stagnation-Point Flow of 
Circular Cylinder in Turbulent Flow 

7. Author(s) 


5. Report Date 


6. Performing Organization Code 


8. Performing Organization Report No. 


Chi R. Wang 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


E— 420 1 

10. Work Unit No. 

505-62-21 

1 1 . Contract or Grant No. 


13. Type o( Report and Period Covered 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


15. Supplementary Notes 


lEMllllHEll 


Memorandum 


14. Sponsoring Agency Code 


Prepared for the 1988 Winter Annual Meeting of the American Society of 
Mechanical Engineers, Chicago, Illinois, November 28 - December 2, 1988. 


16. Abstract 

Boundary layer flow and turbulence transport analyses to study the Influence of 
the free-stream turbulence on the surface heat transfer rate and the skin fric- 
tion around the stagnation point of a circular cylinder in a turbulent flow are 
presented. The analyses are formulated with the turbulent boundary layer equa- 
tions, the Reynolds stress transport equations and the k - e two-equation 
turbulence modeling. The analyses are used to calculate the time-averaged tur- 
bulence double correlations, the mean flow properties, the surface heat transfer 
rate and the skin friction with an Isotropic turbulence in the freestream. The 
analytical results are described and compared with the existing experimental 
measurements. Depending on the free-stream turbulence properties, the turbulence 
kinetic energy can Increase or decrease as the flow moves toward the surface. 
However, the turbulence kinetic energy induces large Reynolds normal stresses at 
the boundary layer edge. The Reynolds normal stresses change the boundary layer 
profiles of the time-averaged double correlations of the velocity and temperature 
fluctuations, the surface heat transfer rate and the skin friction. The free- 
stream turbulence dissipation rate can affect the stagnation-point heat transfer 
rate but the Influence of the free-stream temperature fluctuation on the heat 
transfer rate is insignificant. 

17. Key Words (Suggested by Author(s)) |18. Distribution Statement 


Heat transfer Unclassified - Unlimited 

Boundary layer Subject Category 02 

Turbulence 


19. Security Classlf. (ot this report) 

Unclassified 

20. Security Classlt. (of this page) 

Unclassified 

21. No ol pages 
22 

22. Price" 

A0 2 


NASA FORM 1«se OCT 86 


"For sale by the National Technical Information Service, Springfield, Virginia 22161 
















National Aeronautics and 
Space Administration 

Lewis Research Center 

Cleveland, Ohio 44135 

Official Business 
Penalty tor Private Use $300 


SECOND CLASS MAIL 



LANGLEY BESEARC H CENTER 



1176 01 


327 79 



UESTED 

I 



Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 


NASA 



